Noname manuscript No. 

(will be inserted by the editor) 



Bound states of spin-half particles in a static 
gravitational field close to the black hole field 

A. F. Spencer-Smith • G. H. Gossel • 
J. C. Berengut • V. V. Flambaum 



Received: date / Accepted: date 

Abstract We consider the bound-state energy levels of a spin- 1/2 fcrmion in 
the gravitational field of a near-black hole object. In the limit that the metric 
of the body becomes singular, all binding energies tend to the rest-mass energy 
(i.e. total energy approaches zero). We present calculations of the ground state 
energy for three specific interior metrics (Florides, Soffel and Schwarzschild) 
for which the spectrum collapses and becomes quasi-continuous in the singular 
metric limit. The lack of zero or negative energy states prior to this limit being 
reached prevents particle pair production occurring. Therefore, in contrast to 
the Coulomb case, no pairs are produced in the non-singular static metric. For 
the Florides and Soffel metrics the singularity occurs in the black hole limit, 
while for the Schwarzschild interior metric it corresponds to infinite pressure 
at the centre. The behaviour of the energy level spectrum is discussed in the 
context of the semi-classical approximation and using general properties of the 
metric. 
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1 Introduction 

In this work we consider the problem of a spin-half fermion gravitationally 
bound by a static, massive body described by a non-singular metric. We in- 
vestigate the energy level spectrum of the particle in the limit that the metric 
exhibits a singularity. Previously, we calculated the spectrum of scalar (spin- 
zero) particles bound by such objects pQ and found that it collapses to zero 
energy and becomes quasi-continuous in the singular limit. Importantly, bound 
states with zero energy (where the binding energy is equal in magnitude to 
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the rest mass of the bound particle) were found not to exist when the metric 
is not singular. 

Here we calculate the spin-1/2 fermion energy level spectrum for three 
interior metrics describing the interior of the gravitating (stationary) body: 
the Florides, Schwarzschild, and Soffel metrics, the latter so-named for its 
derivation and subsequent use in [2J. For each metric we calculate the energy 
spectrum of the particle numerically and using approximate analytic expres- 
sions in the near-singular limit of the metric. In the Florides and Soffel space- 
times this limit is characterized by the approach of the Schwarzschild radius 
r s = 2GM/c 2 towards the radius of the body R, where G is the gravitational 
constant and M is the mass of the central body (black hole limit). In the case 
of the Schwarzschild interior metric the singularity occurs for R = 9r s /8 and 
corresponds to a infinite pressure at the origin, rather than the formation of 
an event horizon. 

Our results show that, as in the scalar case, in the limit that the interior 
metric becomes singular, the ground state energy tends to zero; that is, the 
binding energy tends to the rest mass energy of the bound particle. Addition- 
ally, the energy spectrum becomes quasi-continuous in this limit. This means 
that one cannot increase the strength of the central potential (by increasing 
M) to lower a bound state level from the positive continuum to the negative: 
one cannot reach the Dirac Sea to facilitate pair production. This suggests 
that pair-production cannot occur in the field of such a static metric prior to 
the formation of a black hole. 

Additionally, we present a discussion of the semi-classical approach to an- 
alyzing the properties of the level spectrum. It is shown that as the metric 
develops a singularity, all the semi-classical levels will tend to zero. This be- 
haviour can also be understood from classical quantities such as time dilation 
on the interior. 



2 Radial Wave Equations 

One may write a general static, spherically symmetric metric in the form (units 

h=c = i) 

ds 2 = -e^dt 2 + e x ^dr 2 + r 2 dQ 2 . (1) 

The wavefunction for a spin-half particle (Dirac spinor) in a spherically sym- 
metricl field may be separated using 

^=(^^(M)e-- (2) 

with 
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Here k = (— l)' +J+1 / 2 (j + l/2) parametrises the angular momentum for a parti- 
cle with well-defined values of orbital and total angular momentum, j = i±l/2. 
The coupled radial wave equations are then (see, e.g. [2,3 and Appendix A): 



df(r) 
dr 

dg{r) 
dr 



a k 
r 

A K 

r 



+ — f(r) - e» [e *c + m\ g(r) = 0, 



e^—g(r) + e^[e * e — to] f(r) = . 



(4) 
(5) 



These equations can be decoupled to give a single second-order equation 
for the upper component f(r) 
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/(r) = 0. 
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In order to determine the bound states of a spin-1/2 particle in the gravita- 
tional field of a finite-sized object we require a wave function that exists for all 
r. This necessitates an interior as well as an exterior metric to substitute into 
(J6j) and ([2]) . The metric on the exterior of a spherically symmetric non- rotating 
body of mass M is 



(i-r- 



(7) 



where r s = 2GM is the Schwarzschild radius of the gravitating object. This, 
along with a suitable interior metric, allows us to construct interior and exte- 
rior wave equations of the form © which together span all r. These are solved 
both numerically and using approximate analytical formulae. 



3 Method 

3.1 Numeric Approach 

The bound states are found by numerically integrating interior and exterior 
versions of © using Mathematica [7] with the boundary condition for small r 
determined by solving the zero order expansion of the interior wave equation 
around r = 0. The interior wave function is propagated to the boundary where 
continuity of the functions f(R) and f'(R) is imposed. The solution is then 
propagated to large r and the boundary condition f(r —> oo) = is enforced 
by varying e. The s-wave bound states for a given R and to are computed as a 
function of r s . These results are compared with analytic calculations in later 
sections. 



1 



A. F. Spencer-Smith et al. 



3.2 Harmonic Approximation 

In the strong field limit (as the metric becomes singular) the wave function 
of a bound particle is contained entirely on the interior region of the body. 
The assumption of near complete containment on the interior allows one to 
perform all calculations using only the interior wave equation. The validity 
of this assumption is confirmed in subsequent sections. On the interior it is 
convenient to switch to dimensionless variables. We rescale the energy by the 
mass such that s = e/m, and rescale all lengths, including the Compton radius 
of the quantum particle, by the radius of the object R: p = r/R, s = r a /R, 
fx = mR. 

One can re- write ^ in the form of a one dimensional Schrodinger equation 
using the transformation /(/?) = T(p)ip(p) where 

T(p)= (e x /^p{l + ee-/^ ' (8) 

yielding 

-r(p) + v(p)i>( P ) = o (9) 

where V(p) is an effective (energy dependent) potential. 

Analytic calculations of the bound state spectrum in the tight-binding limit 
can be performed by approximating the effective potential on the interior as 
that of a harmonic oscillator. For our interior metrics, one may Taylor 
expand V(p) near the origin and match the resulting expression to that of a 
quantum harmonic oscillator, given by 

V HO (p) = + A^V - (10) 

P 

and use known solutions E n = (n + 3/2)w, where n = 2k + I, k > 0, to 
approximate the energy levels of the system. Here E n = E n (e, s, k, fx) is the 
effective energy of the harmonic oscillator system. 

In the Florides case this approach is valid all the way to the singular 
metric limit. As will be seen in later sections, this is not true for the Soffel and 
Schwarzschild metrics where the harmonic approximation breaks down as the 
metric develops a singularity. 



3.3 Semi-classical approximation 

As the metric develops a singularity, the potential well becomes increasingly 
deep causing the wavefunction to oscillate many times on the interior. In this 
case one may invoke the semi-classical approximation and calculate excited en- 
ergy levels using Bohr-Sommerfield quantization. When treating the problem 
semi-classically, all three singular metric limits may be taken. 
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To do this we investigate Eqn. © for a metric possessing a singularity such 
that e"W ->■ for some value of r. In this limit Eqn. © reduces to 







allowing us to define the approximate semi-classical momentum 

1/2 



P(P) 



(11) 



(12) 



Assuming that the semi-classical approximation is valid, the total phase accu- 
mulated at the boundary will contain the integral 



p(p)dp : 



Pi 



pi 



(A-i/)/2 _ P e 



2s 



dp 



= Xe„-^ 

e 

Setting this phase to kir yields the energy levels as 



irk + v/7r 2 fc 2 + 4(j?XY 
2A> ' 



(13) 



(14) 



This expression tends to zero in the singular limit if Y/X — > 0, which is 
satisfied in this limit as X contains the most singular term. Thus all levels 
tend to zero and the spectrum collapses. 

Additionally, one may arrive at this conclusion by directly considering the 
metric. In the singular limit, the coefficient of dt 2 in the metric tends to zero 
giving rise to (local) time dilation that approaches infinity, with the classical 
period T of a bound particle being similarly affected. Therefore in this limit 
the semi-classical level spacing becomes u> = 2tt/T — > and we arrive at the 
previous result. 



4 Solution of the wave equation in the Florides interior metric 

An interior metric that does not contain the singularity present in the Schwarzschild 
interior solution, discussed in Section ([B]), is that developed by Florides [3]: 



e ^ = (l- S p 2 )- 1 (15) 

which is valid for r ^ R. For r s = R the horizon (coordinate) singularity exists 
at r = r s [5]- The Florides metric corresponds to a spherically symmetric, 
constant density body with vanishing radial stresses. 
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For our purposes this metric is characterized by its effective potential V{p) 
which, near p = (for s — > 1), is given by 

V(p) « ; + M 2 + *« 2 - 7 - T, \T/9 ( 16 ) 



p 2 ^ 4 (l_ s )3/2 



Equating the coefficients of the powers of p between Eqns. (TTT))) and (|16l) . we 
derive two equations for the two unknowns to and e and can therefore solve 
for both. The energy solutions that corresponds to the energy eigenvalues are 



± (1 * )3/ V "f+/3f (17) 



where 



fr=-i»+|)y(»+|) + a r+ 2 «(«-i)-i 

The above expressions are valid under the assumption that the entire wave- 
function is contained within the interior region p < 1 (r < R). This condition 
can be satisfied when /i is large, forcing the particle to be bound near the 
origin. Formally we require (p 2 ) -C 1. The value of (p 2 ) can be obtained using 
the virial theorem for our harmonic potential 

^) = i(n + 3/2) W 
(p 2 ) = « 1. (18) 

/LtCJ 

Expanding a; about l//i in <jXH]) using (IT7|) to lowest order gives u> ~ \/s/2 
where our truncation of the expansion is valid for 



P » + 2P, 



(19) 



recalling that n = 2k + I. Using this value of u> the condition (fT8|) becomes 

M>y^(« + 3/2) (20) 
which automatically satisfies Eqn. (TTHl) . 
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A second test is whether the expansion of Vi(p) as a harmonic potential 
is valid. The ratio of the fourth order to second order term (at the expected 
particle position) is, for large /i, given by 



= — (n + 3/2)V2s 



4/x 



where A n is the coefficient of the p n term in the expansion (|T5|) . Where the 
harmonic oscillator model is appropriate this should be much less than one, 
which leads to 

^> ~(n + 3/2)V2i. (21) 

This condition is automatically satisfied by Eqn. (|20[) . Furthermore, uj > 
is satisfied automatically (when (|2T)|) is satisfied) as it must to be physically 
valid. Finally, we note that in the limit fi > n, we can expand the energy (|17p 
as 

±(l- S ) 3 / 4 (l + (n + 3/2) A /^+o(-^ . (22) 



2^ 2 \fl* 

This is always e < 1 under condition (|2"0"]) (that is, the state is bound). 

Figure ((TJ shows the numeric ground state energy as a function of s com- 
pared with Eqn. <jXTJ) . Figure © compares the analytics and numerics for the 
ground and first three excited states. The figures show that there is good agree- 
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Fig. 1 Ground state energy of spin-1/2 particle in Florides space-time with fj, ■■ 
and k = —1. Circles: numeric; solid line: analytic approximation H17I I. 



10, n = 



ment between our numerical and analytic spectrum. A large mass of [i = 100 
was chosen for Fig. @ in order to satisfy the condition (|2"0")) , which is better 
fulfilled as /i increases, ensuring the validity of the harmonic approximation 
for all n displayed. 
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Fig. 2 Lowest 4 energy levels of spin-1/2 particle in Florides space-time with /i = 100, and 
k = — 1. Circles: numeric; solid line: analytic approximation 1171 . 



5 Solution of the wave equation in the Soffel interior metric 

As an alternative to the approach adopted by Florides, Soffel et al. extended 
the interior Schwarzschild metric past the pressure singularity by analytic 
continuation [5] . We perform our analysis of the bound state spectrum using 
the same procedures and techniques employed in the previous two sections. 
This 'Soffel' metric corresponds to 



Using the methods outlined previously, we find in the Soffel metric the scaled 
energy, e, is 



(1 — s) exp 



. 2(1-*). 



(23) 




(24) 



(l-s) 1 / 2 exp[ 



4(l~s) 



s 



(25) 



e = ± 
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where 



2 + / 2^ («-l)(l-2a) (l-2s)(n + 3/2)< 



4(1 - s) 1 - a 




1 - s 



2 



, o 9/^s o 3 
C = a /r h sk 2 



4(1 - a) 



4 8(1 -a), 
(4 [2s 2 (3k - 4) + k(5 - lis) + 14s] - 26) . (26) 



As in the previous case we examine conditions on [i (taken to be large) 
such that {p 2 } <C 1 and the quadratic term in the expansion of the potential 
dominates (using Eqn. (|2"5l) as e). This yields the condition 



^(n + 3/2) 

which is automatically satisified for s — > 1 (a condition on /i w.r.t. k and 
k can be generated in a similar fashion to Eqn. (|2"0]l of the Florides case, 
but is automatically satisfied by (12"T1) h This is due to the presence of extra 
1/(1 — s) terms in the expansion of the effective potential corresponding to the 
Soffel metric. This also affects the ratio of the fourth order to second order 
coefficients. We find that, unlike in the Florides case, the effective potential 
changes in shape significantly near the origin as s — > 1. Thus, in order to 
ensure that the potential where the particle is bound remains harmonic we 
require a more stringent condition on /i. Specifically, we find that 



(P 2 )A 4 



A 



2 



(n + 3/2)V5 <<: 
2^V2 v / T Tr s 



For comparison with results presented in Fig. Q for n = 4, s = 0.95 this 
corresponds to \l ^> 8. 

As s 1 the energy expression (|2"S")) becomes 



± exp 



4(1 - a) 



which tends to zero in the black hole limit for both positive and negative energy 
states (Eqn. ([28| ensures the second term is <gC 1 in this limit). Figures ([3]) 
and (@} show the ground state and excited state spectra respectively, with 
the numeric and analytic results displayed with each other as before. In the 
singular limit, s — > 1, the spectrum (|25[) is observed to exhibit the same limit 
as the Florides spectrum, (|17[) . The positive energy spectrum collapses to zero 
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Fig. 3 Ground state energy of spin-1/2 particle in Soffel space-time with fi = 10, n = 
and k = —1. Circles: numeric; solid line: analytic approximation (1251) ■ 




Fig. 4 Lowest 4 energy levels of spin-1/2 particle in Soffel space-time with fi 
k = —1. Circles: numeric; solid line: analytic approximation (125 B - 



100. and 



from above as the gravitating body tends towards the black hole limit, and 
the negative energy levels (Dirac sea) also collapse to zero. 

As Eqn. indicates, and as was noted in Sec. (I3.2p . the harmonic ap- 
proximation breaks down in the singular limit for this metric, i.e. only in the 
limit of infinite fi is the harmonic approximation valid. Using the definitions 



Title Suppressed Due to Excessive Length 



11 



in Eqn. (Ti~3l we calculate 



S-J-l 



X = \pK exp 



4(1 - a) 



Y s = 1 (30) 



yielding e — > for all levels as s — > 1. 



6 Solution of the wave equation in the Schwarzschild interior 
metric 

As previously noted, the standard Schwarzschild interior metric cannot be used 
to investigate the limit s — > 1. This is due to a pressure singularity forming at 
the origin for s = 8 /9 where the coefficient of dt 2 is zero. 

The interior solution for a constant density fluid sphere developed by 
Schwarzschild 17 is given by 



* - (I 

e^=(l- S p 2 ) _1 . (31) 



Taking the singular limit (s — > 8/9) of the effective potential yields 
T k(k + 1) / 2 2 sk s(1-2k) 4/z 2 e 2 



p 2 V 2 2(3VT^7-1) (3VT=«- l) 2 

, 2 9s K 2 , s{k-1) 4(2-3yT^)/i 2 e 2 ^ 2 , n , 4 n 



2(371^7 - l) 2 (3VT - 1 - 1) 



(32) 



By equating powers of p between l|10p and (|32p one can solve for e, giving the 
energies as 

(3 % /l — 7-1) 



e = ±^ ■ J -y/a s + 0s, (33) 



where 



S (2-37l~7)(n + 3/2) 2 
" S = 2(3^-1) + ' { ) 



., ,, , S (n + 3/2)2 7 

Sc 



/3s = \/o§c-^+ v ' / ■ (35) 



5=1\u 2 +sk\k-^U, (36) 
7 = ^ + SK ^_9\ (37) 

v s; 2(3yr^~i) 2 
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We can now check the validity of the approximations used. Once again 
invoking the virial theorem and using Eqn. (|18|). the condition (p 2 ) 1 is 
satisfied for 



^ ^ (n + 3/2)V3Vl-3-l s ^8/9 9(71 + 3/2)^8/9-5 



As with Eqn. (|20j) in the Florides condition on /i w.r.t. k and k can be 

generated, but is automatically satisfied by (l38l) . As done previously we require 
the quadratic terms in the expansion to dominate over higher order terms. 
In this case the Schwarzschild metric resembles the Soffel metric considered 
previously, as it too has an additional condition on fi compared to the Florides 
metric. Specifically, we find: 



(P 2 )A 4 



-4, 



(n + 3/2) 



6/^8/9 -s 
(n + 3 
6^/8/9 



« 1, 



,»^#a. 09) 



By taking into account condition (|38p we simplify the energy expression (|33[) 
to 



3^^-l / , (, + 3/2)Vi 

2 V MV^Vl - 7s - 1 J 

which retains the behaviour lim s _j. 8 /ge = for both positive and negative 
energy states. 

Figures ([3]) and ^ show the ground state and excited state spectra re- 
spectively, with the numeric and analytic results displayed with each other as 
before. As in the Florides case, our analytic formula reproduces the bound 
state energies well. In the singular limit, the spectrum is observed to exhibit 
the same qualitative behaviour as that in the Florides case. The positive en- 
ergy spectrum collapses to zero from above, but with one important difference: 
the gravitating body develops a pressure singularity, causing the spectrum to 
collapse, well before the back hole limit. 

As in the Soffel case, the harmonic approximation breaks down when at- 
tempting to model the Schwarzschild case in the s — > 8/9 limit; one must take 
the infinite /i limit to ensure applicability as per Eqn. (|38l) . The corresponding 
calculation in the semi-classical picture is done using 

s^8/9 7T 
J\ — 



v/2(8/9 - 7) ' 

s^8/9 , 

Y « 0.7. (41) 



Substituting this into Eqn. (Tl4"|) yields e — > for all levels as s —> 8/9. 
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Fig. 5 Ground state energy of spin-1/2 particle in Schwarzschild interior space-time with 
fi = 10, n = and n = — 1. The metric becomes singular as s approaches 8/9 causing the 
binding energy to tend towards the rest mass of the fermion. Circles: numeric; solid line: 
analytic approximation (133 D - 
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0.80 0.83 0.85 0.88 
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Fig. 6 Lowest 4 energy levels of spin-1/2 particle in Schwarzschild interior space-time with 
fi = 100, and n = — 1. Circles: numeric; solid line: analytic approximation 1331 . 

7 Discussion 

We present calculations of the energy level spectrum for a massive spin-1/2 
particle bound by a spherical static body. This is done by matching a suitable 
interior metric to the standard Schwarzschild exterior metric and solving the 
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Dirac equation. Numerical computation of the ground state energies of the 
particle shows that the existence of a zero energy bound state may only be 
possible in the limit that the interior metric develops a singularity. In the case 
of the Florides and Soffel metrics this corresponds to the black hole limit; for 
the Schwarzschild interior metric this singularity represents infinite pressure at 
the centre. This result is verified using an approximate analytical calculation. 
For the Florides and Soffel metrics it is seen that the entire bound state spec- 
trum collapses to e = as r s — > R: the energies and the intervals between the 
energy levels are proportional to (1— r s /i?) 3 / 4 and exp[— ■p Tr ^__] respectively. 
If we keep energy fixed the principal quantum number of the level n tends to 
infinity. In the case of the Schwarzschild interior metric the singularity occurs 
for r s — > 8R/9 with the spectrum going as ea \/8/9 — r s /R. 

What is clear is that for any fixed m, n, k, the general form of the equations 
for e n in the Florides and Soffel cases ensures that the bound state energy 
collapses to zero in the limit r s — > R. Additionally, the form of the energy 
spectrum in the Schwarzschild interior case results in the same collapse of 
the spectrum, except now the collapse occurs in the limit r s — > 8R/9. This 
means that for spin-1/2 particles, as in the case of scalar particles, the event 
horizon singularity is not the only singularity that results in collapse of the 
spectrum to a gapless state — a pressure singularity gives the same result. 
Therefore the bound state energy spectrum of an object with a near-singular 
metric reproduce those of an object that is close to the black hole limit. 

The existence of a bound state with zero energy is relevant to the phe- 
nomenon of particle pair production. Due to quantum fluctuations, particle 
anti-particle pairs are produced around a black hole and it is possible that one 
of the pair escapes the gravitational field yielding Hawking radiation [T81IT9] . 

Other systems that give rise to pair production, also known as vacuum 
breakdown, include static and dynamic electromagnetic fields (see, e.g. [20] ) 
and time- varying gravitational fields (see, e.g. [2T] ). In the Coulomb case it 
is clear that one of the pair must be repelled, since the pair have opposite 
charge. In the gravitational case, however, both particles will possess the same 
'gravitational charge' and thus will be attracted by the potential. In the black 
hole case the existence of the event horizon as a barrier facilitates one particle 
escaping to infinity as the other falls to the singularity. 

For gapless energy states to exist in the metrics we have considered, the 
metric must become singular. For pair production with the ejection of one 
particle to infinity it is actually necessary to have a negative energy level 
e < — mc 2 . As a concrete example, compare the gravitational field case with the 
electrostatic potential case where pair creation is possible in a strong Coulomb 
field U(r). (For example, the ground state of an electron orbiting a finite-size 
nucleus reaches the lower continuum when Z > 170 [22]). In the Coulomb case 
we have (e — U(r)) 2 in the wave equation. Therefore, increasing the absolute 
value \U\ of the negative potential U(r) leads to the negative energy e. In 
the gravitational field we have the field dependent term appearing as a 
multiplicative factor of the energy, rather than as a subtractive function, a 
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situation that cannot introduce negative energy bound states. 
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A Derivation of Curved Space Dirac Equation 

Although the following discussion is already well represented in the literature (see, for ex- 
ample \23\), we restate the most important results here as a reference for the reader. 
The Dirac equation on a curved manifold (with ti = c = 1) reads 

(H M V M - m)9 = 0. (42) 

The Dirac matrices 7^ satisfy the Clifford algebra 

T^f+^T" =2gi iv , (43) 

where gV is the reciprocal of the metric and m is the rest mass of the spin-1/2 particle. A 
tilde above a gamma matrix indicates that the matrix is to be taken in curved space-time, 
conversely, the absence of a tilde indicates that the gamma matrix is to be taken in fiat 
Minkowski space-time. The covariant derivative is 

= (9 M - r M )#, (44) 

where the spin-affine connections, i - ^, may be determined from the relation 

MX) = -±y a y pr f> b * (V 8 A - r- A b^) . (45) 

The vierbein of vector fields, are found form the equality 

9y.v = V^b^, (46) 
where r\ a p is the Minkowski metric (with signature ( — M — h) here). The Christoffel symbols 

r^xix) = \g™ Gw,a + g a x,„ - 9 m a, q ) ■ (47) 

This then yields the equations for f(r) and g(r): 

df(r) X k „, » x r " 1 , » , , 

' + c2 -f(r) - C2 |e _ 2 e + m| g(r) = 0, (48) 



dr 

ig(j 

dr r 



dg(r) x k ^ r u 1,, , , 

. e 2 -g(r) + el |e~ 2 e - mj /(r) = 0. (49) 



